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Abstract. Brolin-Lyubich measure \r of a rational cndomorphism R : C — > 
C with degR > 2 is the unique invariant measure of maximal entropy h\ R = 
htop(R) = logd. Its support is the Julia set J(R). We demonstrate that is 
always computable by an algorithm which has access to coefficients of R, even 
when J(R) is not computable. In the case when R is a polynomial, Brolin- 
Lyubich measure coincides with the harmonic measure of the basin of infinity. 
We find a sufficient condition for computability of the harmonic measure of a 
domain, which holds for the basin of infinity of a polynomial mapping, and 
show that computability may fail for a general domain. 



1. Foreword 

This paper continues the line of works [5, 4, 3, 6, 7] of several of the authors on 
algorithmic computability of Julia sets. In this brief introduction we outline our 
results and attempt to give a brief motivation for them. 

Numerical simulation of a chaotic dynamical system: the modern para- 
digm. A dynamical system can be simple, and thus easy to implement numerically. 
Yet its orbits may exhibit a very complex behaviour. The famous paper of Lorenz 
[15], for example, described a rather simple nonlinear system of ordinary differen- 
tial equations x'(t) = F(x) in three dimension which exhibits chaotic dynamics. In 
particular, while the flow of the system $t(a;o) is easy to calculate with an arbitrary 
precision for any initial value xo and any time t. However, any error in estimating 
the initial value xq grows exponentially with t. This renders impractical attempting 
to numerically simulate the behaviour of a trajectory of the system for an extended 
period of time: small computational errors are magnified very rapidly. If we recall 
that the Lorenz system was introduced as a simple model of weather forecasting, one 
understands why predicting weather conditions several days in advance is difficult 
to do with any accuracy. 

On the other hand, there is a great regularity in the global structure of a typical 
trajectory of Lorenz system. As was ultimately shown by Tucker [26], there exists 
a set A C R 3 such that for almost every initial point xq, the limit set of the orbit, 

uj(xo) = A. 

This set is the attractor of the system [25, 19]. Moreover, for any continuous test 
function ?/>, the time average of ip along a typical orbit 
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converges to the integral J ipd/j, with respect to a measure fj, supported on A. 

Thus, both the spatial layout and the statistical properties of a large segment of 
a typical trajectory can be understood, and, indeed, simulated on a computer: even 
mathematicians unfamiliar with dynamics have seen the butterfly-shaped picture of 
Lorcnz attractor A. This example summarizes the modern paradigm of numerical 
study of chaos: while the simulation of an individual orbit for an extended period 
of time does not make a practical sense, one should study the limit set of a typical 
orbit (both as a spatial object and as a statistical distribution). A modern summary 
of this paradigm is found, for example, in the article of J. Palis [21]. 

Julia sets as counterexamples, and the topic of this paper. Julia sets are 
repellers of discrete dynamical systems generated by rational maps R of the Rie- 
mann sphere C of degree d > 2. For all but finitely many points z € C the limit 
of the n-th preimages R~ n (z) coincides with the Julia set J(R). The dynamics of 
R on the set J is chaotic, again rendering numerical simulation of individual orbits 
impractical. Yet Julia sets are among the most drawn mathematical objects, and 
countless programs have been written for visualizing them. 

In spite of this, two of the authors showed in [6] that there exist quadratic 
polynomials f c (z) — z 2 + c with the following paradoxical properties: 

• an iterate f c (z) can be effectively computed with an arbitrary precision; 

• there does not exist an algorithm to visualize J(/ c ) with an arbitrary finite 
precision. 

This phenomenon of non-computability is rather subtle and rare. For a detailed 
exposition, the reader is referred to the monograph [7]. In practical terms it should 
be seen as a tale of caution in applying the above paradigm. 

We cannot accurately simulate the set of limit points of the preimages (f c )~ n (z), 
but what about their statistical distribution? The question makes sense, as for all 
z =/= oo and every continuous test function tjj, the averages 
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where A is the Brolin-Lyubich probability measure [8, 16] supported on the Julia set 
J(/ c ). We can thus ask whether the value of the integral on the right-hand side 
can be algorithmically computed with an arbitrary precision. 

Even if J(/ c ) = Supp(A) is not a computable set, the answer does not a priori 
have to be negative. Informally speaking, a positive answer would imply a dramatic 
difference between the rates of convergence in the following two limits: 
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The main results of the present paper are the following: 

Theorem A. The Brolin-Lyubich measure is always computable. 

The result of Theorem A is uniform, in the sense that there is a single algorithm that 
takes the rational map R as a parameter and computes the corresponding Brolin- 
Lyubich measure. Surprisingly, the proof of Theorem A docs not involve much 
analytic machinery. The result follows from the general computable properties of 
the relevant space of measures. 
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Using the analytic tools given by the work of Drasin and Okuyama [10], or Dinh 
and Sibony [9], we get the following: 

Theorem B. For each rational map R, there is an algorithm A{R) that computes 
the Brolin-Lyubich measure in exponential time. 

The running time of A(R) will be of the form cxp(c(i?) ■ n), where n is the precision 
parameter, and c(R) is a constant that depends only on the map R (but not on n). 
Theorems A and B are not comparable, since Theorem B bounds the growth of the 
computation's running time in terms of the precision parameter, while Theorem A 
gives a single algorithm that works for all rational functions R. 

Lastly, the Brolin-Lyubich measure for a polynomial coincides with the harmonic 
measure of the complement of the filled Julia set. As shown in [6] by two of the 
authors, the filled Julia set of a polynomial is always computable. In view of 
Theorem A, it is natural to ask what property of a computable compact set in the 
plane ensures computability of the harmonic measure of the complement. We show: 

Theorem C. // a closed set K C C is computable and uniformly perfect, and has 
a connected complement, then the harmonic measure of the complement is com- 
putable. 

It is well-known [17] that filled Julia sets are uniformly perfect. Theorem C thus 
implies Theorem A in the polynomial case. Computability of the set K is not 
enough to ensure computability of the harmonic measure: we present a counter- 
example of a computable closed set with a non-computable harmonic measure of 
the complement. 

2. Julia sets of rational mappings 

2.1. Dynamics on the Riemann sphere. We attempt to summarize here for 
convenience of a reader, unfamiliar with Complex Dynamics, the basic facts about 
Julia sets of rational mappings. An excellent book of Milnor [20] presents a detailed 
and self-contained introduction to the subject; proofs of most of the facts we state 
can be found there. 

We first recall that the Riemann sphere C is the Riemann surface with the 
topological type of the 2-sphere, S 2 . Such a complex-analytic manifold can be 
constructed by gluing together two copies of the complex plane C\ — C, Ci = C 
by identifying z G C\ \ {0} with w = 1/z e Ci. This procedure can be loosely 
described as adjoining a point at infinity to the complex plane C\ - we will denote 
co the origin in C2 (so that "00 — 1/0" ) . It is convenient sometimes to visualize C 
as the unit sphere 

S 2 = {x 2 +y 2 + z 2 = 1} cl 3 . 
To this end, consider the stereographic projection from the "north pole" (0, 0, 1) C 
S 2 , which sends S 2 \ {(0, 0, 1)} to the plane z = which we naturally identify with 
Ci = C by z = x + iy. In this model, the north pole becomes the point at infinity. 

The Euclidean metric on R 3 restricted to S 2 is transferred by the stereographic 
projection to the spherical metric on C. This metric is given by 

*■ - (4) 

We will refer to the spherical distance as d^,, as opposed to the usual Euclidean 
distance d. 
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A rational function R(z) = P(z)/Q(z) induces an analytic covering C — > C 
branched at the finitely many critical points £ € C with i?'(C) = 0. The degree d 
of this covering is finite, and coincides with the algebraic degree of R: 

d = max(deg(P),deg(Q)), 

assuming P and Q have no common factors. Every analytic branched covering of 
C of a finite degree is given by a rational function. 

We will consider a rational mapping R of degree degi? = d > 2 (that is, non- 
linear) as a dynamical system on the Ricmann sphere; and denote R n the n-th 
iterate of R. The i?-orbit of a point ( is the sequence {R n (()}^ D =Q . The Julia set is 
defined as the complement of the set where the dynamics is Lyapunov-stable: 

Definition 2.1. Denote F(R) the set of points z E C having an open neighborhood 
U(z) on which the family of iterates R n \u(z) is equicontinuous; that is for every 
e > there exists S > such that if dg.(z, w) < 8 then for every n € N one has 
dfc(R n (z) , R n (w)) < e. The set F(R) is called the Fatou set of R and its complement 
J(R) =C\ F(R) is the Julia set. 

In the case when the rational mapping is a polynomial 

P(z) = a + a x z -\ h a d z d : C ->• C 

an equivalent way of defining the Julia set is as follows. Obviously, there exists 
a neighborhood of oo on C on which the iterates of P uniformly converge to oo. 
Denoting A{o6) the maximal such domain of attraction of oo we have A(oo) C F{R) . 
We then have 

J{P) = dA(oo). 

The bounded set C\A(oo) is called the filled Julia set, and denoted K{P); it consists 
of points whose orbits under P remain bounded: 

K{P) = {zeC\ sup \P n (z)\ < oo}. 

n 

For future reference, let us summarize in a proposition below the main properties 
of Julia sets: 

Proposition 2.1. Let R : C — > C be a rational function. Then the following 
properties hold: 

(a) J(R) is a non-empty compact subset of C which is completely invariant: 
R-^JiR)) = J(R); 

(b) J(R) = J(R n ) for all n e N; 

(c) J(R) has no isolated points; 

(d) if J{R) has non-empty interior, then it is the whole ofC; 

(e) let U C C be any open set with U D J(R) =/= 0. Then there exists n € N such 
that R n (U) D J(R); 

(f) periodic orbits of R are dense in J(R). 

Let us further comment on the last property. For a periodic point zq = R p (zq) of 
period p its multiplier is the quantity A = X(zo) = DR p (zq). We may speak of the 
multiplier of a periodic cycle, as it is the same for all points in the cycle by the Chain 
Rule. In the case when |A| ^ 1, the dynamics in a sufficiently small neighborhood 
of the cycle is governed by the Mean Value Theorem: when |A| < 1, the cycle is 
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attracting (super- attracting if A = 0), if |A| > 1 it is repelling. All repelling periodic 
points are in the Julia set, and all attracting ones are in the Fatou set. 

The situation is much more complicated when |A| = 1; understanding of the local 
dynamics in this case is not yet complete. 

One of the founders of the subject, P. Fatou, has shown that that for a rational 
mapping R with deg R — d > 2 at most finitely many periodic orbits are non- 
repelling. A sharp bound on their number depending on d has been established by 
Shishikura; it is equal to the number of critical points of R counted with multiplicity: 

Fatou-Shishikura Bound. For a rational mapping of degree d the number of the 
non-repelling periodic cycles taken together with the number of cycles of Herman 
rings is at most 2d — 2. For a polynomial of degree d the number of non-repelling 
periodic cycles in C is at most d — 1 . 

Therefore, we may refine the last statement of Proposition 2.1: 
(f) J{R) = {repelling periodic orbits of R}. 
We also note a useful corollary of Proposition 2.1 (e): 

Corollary 2.2. Let w e J(R). Then 

J(R) = |J f- k (w). 

k>0 

To conclude the discussion of the basic properties of Julia sets, let us consider 
the simplest examples of non-linear rational endomorphisms of the Riemann sphere, 
the quadratic polynomials. Every affine conjugacy class of quadratic polynomials 
has a unique representative of the form f c (z) = z 2 + c, the family 

f c (z) = z 2 + c, c e C 

is often referred to as the quadratic family. For a quadratic map the structure of 
the Julia set is governed by the behavior of the orbit of the only finite critical point 
0. In particular, the following dichotomy holds: 

Proposition 2.3. Let K = K(f c ) denote the filled Julia set of f c , and J = J(/ c ) = 
dK. Then: 

• S K implies that K is a connected, compact subset of the plane with 
connected complement; 

• (fi K implies that K = J is a planar Cantor set. 

The Mandelbrot set Ai C C is defined as the set of parameter values c for which 
J(/ c ) is connected. 

2.2. Brolin-Lyubich measure on the Julia set. 

Definition 2.2. Consider a rational map R : C —> C of degree d > 2. We say that 
a probability measure fi on C is balanced (with respect to R) if for every set IcC 
on which R is inject ive we have 

n(R(X))=d-n(X), 



that is, the Jacobian of fi is equal to d. 
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We see that a balanced measure [i is necessarily invariant: as most points in C 
have d preimages under R, 

li(R- 1 (X))=»(X). 

However, a rational map has many invariant probability measures (as a simplistic 
example, for a periodic orbit zq h-> Z\ i-> ■ ■ ■ z v -\ H> zq define n = |X^z<)- On 
the other hand there is exactly one balanced measure for R: the Brolin-Lyubich 
measure A. Constructed by Brolin [8] for polynomials, and later by Lyubich [16] for 
a general rational function, it is supported on the Julia set J(R). Lyubich showed 
that for all but finitely many points zsC the weak limit 

lim i E 5 ^ = x - 

n— >oo a * — * 

In general, given a transformation T of a compact space X, denote by h top (T) 
and hn(T) the topological and measure-theoretic entropies, respectively. The well 
known Variational Principle, tell us that: 

hop(T) = sup h^T), 

[iEMt 

where Mt denotes the set of T-invariant measures. A measure n is called a measure 
of maximal entropy if h^(T) = h top (T). 

Lyubich showed that A is the unique measure of maximal entropy of R: 

Theorem 2.4 ([10]). The measure A is the unique measure on C for which the 
metric entropy h\(R) coincides with the topological entropy of R: 

h x (R) = h top (R) = logd. 

Note that for any invariant measure \x we have 

J Jac M /cfyi < d, 

therefore a measure of maximal entropy is necessarily balanced. 

2.3. Harmonic measure in polynomial dynamics. A detailed discussion of 
harmonic measure can be found in [12]. Here we briefly recall some of the relevant 
facts. 

Let G be a simply-connected domain in C whose complement K contains at least 
two points, and g £ G. The harmonic measure uJc,g is defined on the boundary dG. 
For a set E C dG it is equal to the probability that a Brownian path originating 
at g will first hit dG within the set E. 

To define the harmonic measure for a non simply-connected domain G = C \ K 
we have to require that a Brownian path originating in G will hit dG almost surely, 
a condition which is satisfied automatically for a simply-connected domain. A 
quantitative measure of a likelyhood that such a set will be hit by a Brownian path 
is defined as follows. Consider K <s C, and let B t be a Brownian path which is 
started uniformly at a circle {|z| = R} which surrounds K. Denote r the first 
moment when B T £ K . The logarithmic capacity of K is 



CapOK) =exp(E(log|B T |)). 
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By way of an example, consider a connected and locally-connected compact set 
K C C. In this case, dG is a continuous image of the unit circle. In fact, consider 
the unique conformal Ricmann mapping 

if) : G = C \ K -> C \ D R (0) = (Ar(0)) c , with ifj(oo) = oo and ^'(oo) = 1. 

The quantity r{G, oo) = 1/R is the conformal radius of G about oo. 

By a classical theorem of Caratheodory, ip -1 extends continuously to map G — > 
C \ Dr(0). By symmetry considerations, the harmonic measure Wm R (o)) c ,oo coin- 
cides with the Lcbcsgue measure ^ on the circle dDn(0) = {\z\ = R}. Conformal 
invariancc of Brownian motion implies that uj£j ;00 is obtained by pushing forward 

M by V^lazMO); and that 

Cap(G) = l/r(G,oo). 

Consider a polynomial P : C — > C with degP > 2. The capacity of the filled 
Julia set K(P) is equal to one. This follows from a classical result of Bottcher when 
K(P) is connected (see [8] for the general case). Brolin [8] was the first to show that 
the balanced measure A of P coincides with the harmonic measure ^c\K(P) oo- 

3. COMPUTABILITY 

3.1. Algorithms and computable functions on integers. The notion of an 
algorithm was formalized in the 30's, independently by Post, Markov, Church, and, 
most famously, Turing. Each of them proposed a model of computation which de- 
termines a set of integer functions that can be computed by some mechanical or 
algorithmic procedure. Later on, all these models were shown to be equivalent, 
so that they define the same class of integer functions, which are now called com- 
putable (or recursive) functions. It is standard in Computer Science to formalize 
an algorithm as a Turing Machine [27] . We will not define it here, and instead will 
refer an interested reader to any standard introductory textbook in the subject. It 
is more intuitively familiar, and provably equivalent, to think of an algorithm as a 
program written in any standard programming language. 

In any programming language there is only a countable number of possible algo- 
rithms. Fixing the language, we can enumerate them all (for instance, lexicograph- 
ically). Given such an ordered list (A n )%f =1 of all algorithms, the index n is usually 
called the Godel number of the algorithm A n . 

We will call a function / : N — > N computable (or recursive), if there exists an 
algorithm A which, upon input n, outputs f(n). Computable functions of several 
integer variables are defined in the same way. 

A function / : W — > N, which is defined on a subset W C N, is called partial 
recursive if there exists an algorithm A which outputs f(n) on input n € W, and 
runs forever if the input n ^ W. 

3.2. Time complexity of a problem. For an algorithm A with input w the 
running time is the number of steps A makes before terminating with an output. 
The size of an input w is the number of dyadic bits required to specify w. Thus 
for w € N, the size of w is the integer part of log 2 w. The running time of A is the 
function 

T A : N -> N 

such that 

T/i(n) = maxjthc running time of A(w) for inputs w of size n}. 
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In other words, TU(n) is the worst case running time for inputs of size n. For a 
computable function / : N — > {0, 1} the time complexity of / is said to have an 
upper bound T(n) if there exists an algorithm A with running time bounded by 
T(n) that computes /. 

3.3. Computable and semi-computable sets of naturals numbers. A set 

E C N is said to be computable if its characteristic function \e ■ N — > {0, 1} is 
computable. That is, if there is an algorithm A : N — > {0, 1} that, upon input n, 
halts and outputs 1 if n G E or if n ^ E. Such an algorithm allows to decide 
whether or not a number n is an element of E. Computable sets are also called 
recursive or decidable. 

Since there are only countably many algorithms, there exist only countably many 
computable subsets of N. A well known "explicit" example of a non computable 
set is given by the Halting set 

H := {i such that Ai halts}. 

Turing [27] has shown that there is no algorithmic procedure to decide, for any 
ieN, whether or not the algorithm with Godcl number i, Ai, will eventually halt. 

On the other hand, it is easy to describe an algorithmic procedure which, on 
input i, will halt if i G H, and will run forever if i ^ H . Such a procedure can 
informally be described as follows: on input i emulate the algorithm At; if Ai halts 
then halt. 

In general, we will say that a set E C N is lower-computable (or semi- decidable, 
or recursively enumerable) if there exists an algorithm Ae which on an input n 
halts if n G E, and never halts otherwise. Thus, the algorithm Ae can verify the 
inclusion n G E, but not the inclusion n G E c . We say that Ae semi-decides 
n G E (or semi-decides E). The complement of a lower-computable set is called 
upper- computable. 

The following is an easy excercisc: 

Proposition 3.1. A set is computable if and only if it is simultaneously upper- 
and lower-computable. 

3.4. Computability over the reals. Strictly speaking, algorithms only work on 
natural numbers, but this can be easily extended to the objects of any countable 
set once a bijection with integers has been established. The operative power of 
an algorithm on the objects of such a numbered set obviously depends on what 
can be algorithmically recovered from their numbers. For example, the set Q of 
rational numbers can be injectivcly numbered Q = {qo, q\, . . .} in an effective way: 
the number i of a rational a/b can be computed from a and b, and vice versa. 
The abilities of algorithms on integers are then transferred to the rationals. For 
instance, algorithms can perform algebraic operations and decide whether or not 
qi > qj (in the sense that the set {(i,j) ■ qi > qj} is decidable). 

Extending algorithmic notions to functions of real numbers was pioneered by 
Banach and Mazur [1, 18], and is now known under the name of Computable Anal- 
ysis. Let us begin by giving the definition of a computable real number, going back 
to the seminal paper of Turing [27]. 

Definition 3.1. A real number x is called 

• computable if there is a computable function / : N — > Q such that 

|/(n)-x|<2-"; 
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• lower- computable if there is a computable function / : N — > Q such that 

fin) / x; 

• upper- computable if there is a computable function / : N — > Q such that 

f(n) \ x. 

Algebraic numbers or the familiar constants such as 7r, e, or the Feigembaum 
constant are all computable. However, the set of all computable numbers Re is 
necessarily countable, as there are only countably many computable functions. 

We also remark that if x is lower-computable then there is an algorithm to semi- 
decide the set {qi < x}: just compute f(n) for each n and halt if qi < f{n). In 
other words, the set {q E Q : q < x} is lower-computable. The converse is also 
obviously true: 

Proposition 3.2. If E a Q is lower- computable and x = supE 1 < oo, then x is 
lower-computable. 

In the same way as there exist lower-computable sets which are not computable, 
there exists lower-computable numbers which are not computable. The usual con- 
struction is as follows: let {a{)i be an algorithmic enumeration (without repetitions) 
of a lower-computable set A which is not computable. For instance, we can take 

A = {i 6 N such that Ai halts}. 

Define 

</, f> ° ' 

i=0 

Clearly, (q n ) n is a computable non-decreasing sequence of rational numbers. Be- 
ing bounded by 1, it converges. The limit, say x, is then a lower-computable num- 
ber. It x were computable, it would be possible to compute the binary expansion 
of x which, in turn, would allow to decide the set A. 

We also note: 

Proposition 3.3. A real number x is computable if and only if it is simultaneously 
lower- and upper- computable. 

Proof. Let us assume that x is both lower- and upper-computable. Thus there exist 
algorithms A\ and A2 which compute sequences of rationals aj and bj respectively 
with 

aj /*■ x and bj \ x. 

Consider the algorithm A which on the input n emulates A±, A2 to find the first 
k(n) such that \a^ n ) — 6fe( ra )| < 2~™, and then outputs afc(„j. Then f(n) = a^ n ^ is 
a computable function such that \ f(n) — x\ < 2~ n and hence x £ Rc- 

The other direction is trivial. □ 

3.5. Uniform computability. In this paper we will use algorithms to define com- 
putability notions on more general objects. Depending on the context, these objects 
will take particulars names (computable, lower-computable, etc..) but the defini- 
tion will always follow the scheme: 

an object x is computable if there exists an algorithm A satisfying the property 
P(A, x). 
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For example, a real number x is computable if there exists an algorithm A which 
computes a function / : N — > Q satisfying \f(n) — x\ < 2~ n for all n. Each time 
such definition is made, a uniform version will be implicitly defined: 

the objects {x 7 } 7e r are computable uniformly on a countable set T if there exists 
an algorithm A with an input 7 £ T, such that for all 7 £ T, A-y := -4(7, •) satisfies 
the property P(A 1 ,x 1 ). 

In our example, a sequence of reals (xi)i is computable uniformly in i if there exists 
A with two natural inputs i and n which computes a function f(i, n) : N x N — > Q 
such that for all i € N, the values of the function fi(-) :~ f(i,-) satisfy 

|/i(n) -Xi\< 2~ n for all n £ N. 

3.6. Computable metric spaces. The above definitions equip the real numbers 
with a computability structure. This can be extended to virtually any separable 
metric space, making them computable metric spaces. We now give a short intro- 
duction. For more details, see [28]. 

Definition 3.2. A computable metric space is a triple {X, d, S) where: 

(1) (X, d) is a separable metric space, 

(2) S = {si : i £ N} is a dense sequence of points in X, 

(3) d(si,Sj) are computable real numbers, uniformly in (i,j). 

The points in S are called ideal. 

Example 3.1. A basic example is to take the space X = W l with the usual notion 
of Euclidean distance d(-, ■), and to let the set iS consist of points x = (xi, ■ ■ ■ , x n ) 
with rational coordinates. In what follows, we will implicitly make these choices of 
S and •) when discussing computability in R". 

Definition 3.3. A point x is computable if there is a computable function / : N — > 
N such that 

\ s f(n) — x\ < 2~ n for all n. 

If x 6 X and r > 0, the metric ball B(x, r) is defined as 

B(x,r) = {y e X : d(x,y) < r}. 

Since the set B := {B(s, q) : s € S, q £ Q, q > 0} of ideal balls is countable, we can 
fix an enumeration B = {Bi : i £ N}. 

Proposition 3.4. A point x is computable if and only if the relation x £ Bi is 
semi-decidable, uniformly in i. 

Proof. Assume first that x is computable. We have to show that there is an algo- 
rithm A which inputs a natural number i and halts if and only if x £ Bi . Since x is 
computable, for any n we can produce an ideal point s n satisfying |s„ — x\ < 2~ n . 
The algorithm A work as follows: upon input i, it computes the center and radius 
of Bi, say s and r. It then searches for n £ N such that 

d{s ni s) + 2- n <r. 

Evidently, the above inequality will hold for some n if and only if x £ Bi. 

Conversely, assume that the relation x £ Bi, s scmi-dccidablc uniformly in i. To 
produce an ideal point s n satisfying |s„ — x\ < 2~ n , we only need to enumerate all 
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ideal balls of radius 2 71+1 until one containing x is found. We can take s n to be 
the center of this ball. □ 

Definition 3.4. An open set U is called lower- computable if there is a computable 
function / : N — > N such that 

U=\jB f(n) . 

Example 3.2. Let e > be a lower-computable real. Then the ball B(0, e) is 
a lower-computable open set. Indeed: B(so,e) = \J n B(0,q n ), where (<?„)„ is the 
computable sequence converging to e from below. 

It is not difficult to see that finite intersections or infinite unions of (uniformly) 
lower-computable open sets are again lower computable. As in Proposition (3.4), 
one can show that the relation x G U is semi-decidable for a computable point x 
and an open lower-computable set. 

We will now introduce computable functions. Let X' be another computable 
metric space with idea balls B' = {B^}. 

Definition 3.5. A function / : X — > X' is computable if the sets f~ 1 (B' i ) are 
lower-computable open, uniformly in i. 

An immediate corollary of the definition is: 

Proposition 3.5. Every computable function is continuous. 

The above definition of a computable function is concise, yet not very transpar- 
ent. To give its e — 6 version, we need another concept. We say that a function 
4> : N — > N is an oracle for x G X if 

d(s^ (m) ,x) < 2-' m . 

An algorithm may query an oracle by reading the values of the function <f> for an 
arbitrary n G N. We have the following: 

Proposition 3.6. A function f : X — > X' is computable if and only if there exists 
an algorithm A with an oracle for x G X and an input n G N which outputs s' n G S' 
such that d(s' n , f(x)) < 2~". 

In other words, given an arbitrarily good approximation of the input of / it is 
possible to constructively approximate the value of / with any desired precision. 

3.6.1. Computability of closed sets. Having dfined lower-computable open sets, we 
naturally proceed to the following definition. 

Definition 3.6. A closed set K is upper- computable if its complement is lower- 
computable. 

Let us look at two examples. Firstly, 

Example 3.3. A closed ideal ball c\(B(s,q)) is clearly upper-computable. To see 
this, observe that a point s n belongs to X \ c\(B(s,q)) if and only if d(s n ,s) > 
q. Since this last relation is semi-decidable, we can enumerate such ideal points. 
Moreover, for each of them we can also find q n satisfying < q n < d(s ni s) — q, so 
that B(s n ,q n ) cX\cl(B(s,q)). 

Our second example is more interesting: 
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Example 3.4. Let P : C — > C be a computable polynomial of degree degP > 2. 
Then the filled Julia set K(P) is upper-computable. 

Proof. Indeed, let M € Q be such that If (P) C B(0, M). Enumerate the points in 
R 2 with rational coordinates S = {s n = (a n , &«)}, and set £„ = a n + ib n . For every 
point C, n € C \ K(P) we can identify an iterate P z (Cn) ^ 5(0, M). Moreover, for 
such a point we can find e„ € Q such that 

P*(P(C„,e„))nP(0,M) = 0. 

We can thus algorithmically enumerate a sequence of open ideal balls which ex- 
hausts R 2 \ K{P). □ 

Definition 3.7. A closed set K is lower- computable if the relation K Bi 7^ is 
scmi-decidable, uniformly in i. 

In other words, a closed set K is lower-computable if there exists an algorithm A 
which enumerates all ideal balls which have non-empty intersection with K . 

To see that this definition is a natural extension of lower computability of open 
sets, we note: 

Example 3.5. 

(1) The closure of an ideal ball c\(B(s, q)) is lower-computable. Indeed, B(si,qijD 
c\(B(s, q)) ^ if and only if d(s, s n ) < q + q n . 

(2) More generally, the closure cl(t/) of any open lower-computable set U is 
lower-computable since Bi Dc\(U) 7^ if and only if there exists s £ Bi n U. 

The following is a useful characterization of lower-computable sets: 

Proposition 3.7. A closed set K is lower-computable if and only if there exists a 
sequence of uniformly computable points Xi € K which is dense in K . 

Proof. Observe that, given some ideal ball B = B(s, q) intersecting K, the relations 
cl(Pj) C B, qi< 2- k and B t n K ^ are all semi-decidable and then we can find 
an exponentially decreasing sequence of ideal balls (Pfe) intersecting K. Hence 
{x} = HfcPfc is a computable point lying in B n K. 

The other direction is obvious. □ 

Example 3.6. Let R be a computable rational map of degree degP > 2. Then 
the Julia set J{R) is lower-computable. 

Sketch of proof. We will use Corollary 2.2. Periodic points of R are computable (by 
any standard root-finding algorithm) and so are their multipliers. We can semi- 
decide whether a periodic point is repelling (if the multiplier is greater than 1 we will 
be able to establish this with a certainty by computing the point and its multiplier 
precisely enough). Therefore, the repelling periodic points of R are computable. 
Let w be any such point. The points in Uk>of~ k (w) are uniformly computable, 
and dense in J(R). By Proposition 3.7, J(R) is a closed lower-computable set. □ 

Definition 3.8. A closed set is computable if it is lower and upper computable. 

Putting together Examples 3.4 and 3.6, we obtain the following theorem of [4]: 

Example 3.7. Let P be a computable polynomial with degP > 2, and suppose 
that K(P) has empty interior, that is, K(P) = J(P). Then K(P) is a computable 
set. 
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(c) 



FIGURE 1. Computing the Julia set of f(z) = z 2 + i (J(f) = 
K(f)). (a) approximating the filled Julia set from above: the first 
15 preimages of a large disk D = B(0, R) 2> K(f)\ (b) approximat- 
ing the Julia set from below: Uo<fc<i2/~ fc (/3) where j3 is a repelling 
fixed point in J(/); (c) a good-quality picture of </(/). 

Here is an alternative way to define a computable set. Recall that Hausdorff 
distance between two compact sets Ki, K 2 is 

dist^i, i<r 2 ) = inf{^! C U,(K 2 ) and K 2 C U e (K\)}, 

e 

where U C (K) = [J zeK B(z, e) stands for an e-neighborhood of a set. The set of all 
compact subsets of M equipped with Hausdorff distance is a metric space which 
we will denote by Comp(A/). If M is a computable metric space, then Comp(M) 
inherits this property; the ideal points in Comp(M) arc finite unions of closed ideal 
balls in M. We then have the following: 

Proposition 3.8. A set K (s M is computable if and only if it is a computable 
point in Comp(M) . 

Equivalenly, K is computable if there exists an algorithm A with a single natural 
input n, which outputs a finite collection of closed ideal balls cl B% , . . . , cl Bi n such 
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that 

distjydJclB^iT) < 2"™. 

3.7. Computable probability measures. Let M.(X) denote the set of Borel 
probability measures over a metric space X. We recall the notion of weak conver- 
gence of measures: 

Definition 3.9. A sequence of measures \i n £ M(X) is said to be weakly convergent 
to fj, £ M(X) if / fdfi n -> / fdfx for each / £ C {X). 

Any smaller family of functions characterizing the weak convergence is called 
sufficient. It is well-known, that when X is a compact separable and complete 
metric space, then so is M.(X). 

Weak convergence on M(X) is compatible with the notion of Wasserstein-Kantorovich 
distance, defined by: 

Wi{n,v) = sup / fdfi- / fdv 

/Gl-Lip(X) J J 

where l-Lip(A) is the space of Lipschitz functions on X , having Lipschitz constant 
less than one. 

The following result (see [13]) says that, when X is a computable metric space, 
M.(X) inherits its computability structure. 

Proposition 3.9. Let T> be the set of finite convex rational combinations of Dirac 
measures supported on ideal points of X. Then the triple (M. (X) , Wi , T>) is a 
computable metric space. 

Definition 3.1. A computable measure is a computable point in (A4(X), W\, Tf). 
That is, it is a measure which can be algorithmically approximated in the weak 
sense by discrete measures with any given precision. 

As examples of computable measures, consider the Lebesgue measure in R™, or any 
smooth measure in R™ with a computable density function. 

The following proposition (see [13]) gives a useful characterization of the com- 
putability of the measure. 

Proposition 3.10. Let /i be a Borel probability measure. The following statements 
are equivalent: 

(1) fj. is computable, 

(2) //(-B^ U . . . U Bi n ) are lower-computable, uniformly in i\, . . . , i n , 

(3) for any uniformly computable sequence of functions (fi)i, the integral J fidfx 
is computable uniformly in i. 

We will also need the following fact (see [24]), : 

Proposition 3.11. If (fi)i is a uniformly computable sequence of functions, then 
the integral operators 

L t : M(X) -> R defined by Li(ff) := J fidfi, 

are uniformly computable. 
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To illustrate the concepts we have introduced, we end this section by constructing an 
example of a computable set E C / = [0, 1] such that Lebesgue measure restricted 
to it is not computable. Indeed, any non-atomic probability measure assigning 
positive measure to intervals in E, will not be computable. 

Example 3.8. Let Ai be a Godel ordering of algorithms. Set oq = and ai = 
a^_i + 2~ l , Define the set Si to be: 

J (at, aj+i) \ {ai + n2~ j " 1 : n = 1, 2 J — 1}, Ai halts in j steps 
0, docs not halt. 

Our set is defined by 

e = i\ua. 

Clearly, Uj5j is lower-computable open and thus E is upper computable. 

Let us prove that E is also lower-computable by producing a dense computable 
sequence of points in E. To this end, we run an algorithm which at step j simulates 
all algorithms Ai, i < j for the first j steps (or until they halt). For every i such 
that Ai does not halt in fewer than j steps it then outputs the set 

{a i + n2~ j ' i :n = 1, 2 j - 1}. 

We denote Ej the union of the sets output by the algorithm at step j. 
It is clear that 

£7 = ol(U^)- 
Thus E is lower-computable, and hence, computable. 

Suppose [i is a non atomic probability measure on E assigning positive mass to 
every interval in E. Then, for each i, (J,(ai, a,+i) > if and only if Ai does not halt. 

Let us assume that /i is a computable measure on E. Then, by Theorem 3.10, 
the relation > is semi-decidable for any rational interval /. Hence the 

Halting set is upper-computable. Since it is also lower-computable, the Halting set 
is computable. We have thus arrived at a contradiction with the undccidability of 
the Halting problem. 



4. COMPUTABILITY OF BROLIN-LYUBICH MEASURE 

4.1. Some preliminaries. In what follows we will require the following facts. The 
first theorem is classical, see e.g. [20]. 

Koebe 1/4 Theorem. // / is a univalent function on a disk D = B(zo,r) C C, 
then 

dist( f(z ),d(f(D)) > \\f{z )\r 

Considerations of compactness (see [11]) imply that there is an algorithmic pro- 
cedure to semi-decide whether a given lower-computable open set of probability 
measures on C contains the whole M.(C). It will be convenient for us to use a 
uniform version of this statement: 

Proposition 4.1. Let (Ui)i be a sequence of uniformly lower- computable open sub- 
sets of M(C). Then the relation M(<C) C Ui is semi-decidable, uniformly in i. 
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Sketch of proof. It is enough that for any given finite list of ideal balls {-Bfci}™i, 
we can semi-decide the relation 

M(C) C |J B ki . 

If this last relation holds, then the union on the right must contain the elements of 
any 2 _n -net of AA(C), provided that 2~™ is less than (half of) the Lebesgue number 
of the covering {-Ba,-}™ !■ Such a net can be computed from a net of C and a net 
of [0,1]. 

□ 

4.2. Proof of Theorem A. Consider a rational map R(z) — P(z)/Q(z) of degree 
d. The coefficients of P and Q form two (d + l)-tuples of complex numbers, and 
we can thus specify R by a (2d + 2)-tuple of coefficients, or a point in <C 2d+2 . It is 
clear that 

Proposition 4.2. R : C — >■ C is a computable function if and only if there exists a 
computable point in C 2d+2 which specifies R. 

Let us now formulate a precise version of the Theorem A: 

Theorem 4.3. For a rational map R denote by \r its Brolin-Lyubich measure. 
Then the functional 

T : C 2d+2 -> M(C) 
R^X R 

is computable. 

Remark 4.4. In other words, there exists an algorithm A with an oracle for v £ 
C 2d+2 and a single natural input n which outputs a measure /i € D which has the 
following property. If R is the rational map with coefficients v then 

W 1 (n,X R )<2- n . 

Of course, if R is computable, then the oracle can be replaced with an algorithm 
computing the coefficients of R. 

Proof of the Theorem A. Let R be a rational map of degree d and <f> be an oracle 
for the coefficients of R. Given n € N, we will show how to compute an ideal ball 
B n C jM(C) with radius 2~™ containing A^. 

Let U C M. (C) be the set of probability measures which are not invariant with 
respect to R, and let V C Ai(C) be the set of probability measures which are not 
balanced. In the following, we show that, using the oracle <j>, both U and V are 
lower-computable open sets. 

Let us introduce a certain fixed, enumerated sequence of Lipschitz computable 
functions which we will use as test functions. Let Ho be the set of functions of the 
form: 

^ ir) e = |l-|d(x,s)-r| + /e|+ (4.1) 

where s is a rational point in C, r, e € Q and \a\ + = max{a, 0}. These are uniformly 
computable Lipschitz functions equal to 1 in the ball B(s, r), to outside £?(s, r + e) 
and with intermediate values in between. 
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Let 

H = {<pi,<pa,...,} (4.2) 

be the smallest set of functions containing 7io and the constant 1, and closed under 
max, min and finite rational linear combinations. Clearly, we have: 

Proposition 4.5. % is a sufficient family of uniformly computable functions. 

Moreover, the functions inH are of the form ip n — c n + g n where c„ is a constant 
computable from n, and g n has a bounded support, and from n one can compute a 
bound for its diameter. 

Lemma 4.6. The set 

U := {/x € Ai(C) : fi is not invariant with respect to R} 
is a lower- computable open set. 

Proof. We show that U is lower-computable open by exhibiting an algorithm to 
semi-decide whether a probability measure /i belongs to U. By Proposition 3.11 
the numbers 



<Pi d/i and / ipi o Rd/i 
are uniformly computable. If [i is not invariant, then there exist j such that 

gid[i^ J 9i° Rdfi 

and such a j can be found. □ 
Lemma 4.7. The set 

V := {(j, € M(C) : [i is not balanced with respect to R} 
is a lower- computable open set. 

Proof. To semi-decide whether a measure [i is not balanced, we start by enumer- 
ating all the ideal points Zi in C which are not critical for R. For each Zi, we can 
compute an ideal ball Bi = B(zi,r) such that J?|s j is injective. Denote Criti? is 
the set of critical values of R. Compute a rational number q such that 

< q < dist(i?(zj),Crit fl ). 

The function R has a conformal inverse branch RT~ on B(R(zi), q). Compute any 
rational number r such that 

< r < 



4R'(zO 
By Koebe 1/4 Theorem, 

Bi = B(zi,r) C R^ 1 (B(R(zi), q)), 

so that R is conformal on Bi. 

Now, for each Bi, let (tfj)j be the list of test functions supported on Bi. If is 
not balanced, then there exists i such that 

n(R(Bi))^d-fi(Bi) 

which means that there exists / such that 



j if] o RT^dfX + d- J 



ipidfi. 
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By Proposition 3.11, the numbers f <Pj ° R t 1 d/j, and d J <Pjd[j, are uniformly com- 
putable, and thus I can be found. □ 

It follows that the open set U — U U V of measures which are either not in- 
variant or not balanced is lower-computable with an oracle <fi. Its complement is 
the singleton {A#}. To compute A# with precision 2~™, enumerate all the ideal 
balls B n C A4(C) of radius 2~™ and semi-decide the relation {Xr} C B n . This is 
possible because 

{Xr} c B n <^> M(C) cMU B n , 
and the last relation is semi-decidable by Proposition 4.1. 

□ 

4.3. A comparison of rates of convergence. In [6] two of the authors have 
shown: 

Theorem 4.8. There exists a computable quadratic polynomial f c (z) = z 2 + c 
whose Julia set J c is not computable. 

Together with Theorem A this statement has the following amusing consequence: 

Theorem 4.9 (Incommensurability of rates of convergence). For a polyno- 
mial f c and z € C denote 

fa(n) = distff((/ c )~"(z), J c ) and fa(n) = Wi I — S w , A c 

where A c is the Brolin-Lyubich measure of f c . Even though both (f>i and fa converge 
to as n — > oo, there exists a parameter c such that there does not exist any 
computable function F:l4l such that F(0) = and 

fa(n)<Fofa(n). 

Proof. If c is computable, Theorem A implies the computability of A c . Hence, fa(n) 
is a computable function. On the other hand, if there exists a computable bound 
fa{n) — > 0, then J c is a computable set. Therefore, such a bound cannot exist 

n — ^ oo 

for a parameter c as in Theorem 4.8. 

□ 

As an illustration, consider Figure 2. The Julia set of a quadratic polynomial is 
rendered in gray. This particular polynomial can be written in the form Pg(z) — 
z 2 + cxp(2ni9)z for 9 = (\/5 + l)/2. The preimage {P g )~ 12 {z) (highlighted in 
black) for a point z S J(Pe) gives an excellent approximation of A, but a very poor 
approximation of the whole Julia set. 

4.4. Proof of Theorem B. For a given point z € C, set 

^ z > m = ^7 ^ 

The following result is due to Drasin and Okuyama [10], and, in more generality, 
to Dinh and Sibony [9]. 
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Theorem 4.10 ([10, 9]). For each R there are constants a = a(R) and A = A(R) 
such that for every point z£C, except at most two, and for every / el - Lip the 
following holds: 

J fdX M , m -J fdX < Aa~ m . 
Note that A and a are independent of n,z and f. 

We have then that, taking m = Cn for some constant C = C(A, a), 

Wi(A,, m ,A)<2- n . 

So that in order to compute a 2~™-approximation of Aij, it is enough to compute 
approximations of the d m pre- images of z by R m . Since each pre-image can be 
computed in time polynomial in d m , the entire computation can be achieved in 
time <3(2 cn ) for a c = c(R). 

4.5. A counter-example. In view of the above results, it is natural to ask whether 
a measure of maximal entropy of a computable dynamical system is always com- 
putable. The example below will show that this need not be the case. We will 
construct a map 

T : S 1 x S 1 -> S 1 x S 1 

with the following properties: 

(a) T is a computable function; 

(b) T has a measure of maximal entropy; 

(c) every measure of maximal entropy of T is non-computable. 
We first recall a construction [11]: 

Proposition 4.11 ([11]). There exists a computable transformation T\ : S 1 — > S 1 
for which every invariant measure is non- computable. 

To prove this, we need the following facts: 



20 ILIA BINDER, MARK BRAVERMAN, CRISTOBAL ROJAS, MICHAEL YAMPOLSKY 



Proposition 4.12. Let /i be a computable Borel measure on a computable metric 
space X. Then the support of [i contains a computable point x € X. 

Sketch of proof. We outline the proof here and leave the details to the reader. First, 
for each ideal ball B = B(x, r) set 

i'B = 4>x.r/2.r/2 

as in (4.1). An exhaustive search can be used to find a sequence of ideal balls 
B' 1 = B(xi,ri) with the following properties: 

• n < 2-<; 

• B(xi,ri) C B(xi- 1} 2ri_i); 

• J i^Bid^ > 0. 

The algorithm can then be used to compute x = limx, G supp(/x). 

□ 

Proposition 4.13. There exists a lower- computable open set V C (0,1) such that 
(0, 1) \ V 7^ and V contains all computable numbers in (0, 1). 

Sketch of proof. Consider an algorithm A which at step m emulates the first m 
algorithms A%{i), i < m with respect to the Godcl ordering for m steps. That is, 
the i-th algorithm in the ordering is given the number i as the input parameter. 
From time to time, an emulated algorithm Ai(i) may output a rational number Xi 
in (0, 1). Our algorithm A will output an interval 

£, = (^-3-72,^ + 3-72)0(0,1) 

for each term in this sequence. The union V = ULi is a lower-computable set. It is 
easy to see from the definition of a computable real that V D Rc PI (0, 1). If a; € M.c 
then there is a machine A n {j) that on input j outputs a S^^-approximation of x. 
Thus the execution of A n (n) will halt with an output q such that \x — q\ < 3 _n /4, 
and x will be included in V. On the other hand, the Lebesgue measure of V is 
bounded by 1/2, and thus does not cover all of [0, 1]. □ 

Sketch of proof of Proposition 4-11- By Proposition 4.13, there exists an open lower- 
computable set V C [0, 1] such that the complement K = [0, 1] \ V contains no 
computable points. Since V is lower-computablle, there are computable sequences 
{cti, bi}i>i such that < a% < b{ < 1 and V = Ui( a i7 




FIGURE 3. Left: a map right: the map T\ 
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Let us define non-decreasing, uniformly computable functions /» : [0, 1] 
such that 

fi{x) > x if x <E (oj, hi) and fi(x) = x otherwise. 
For instance, we can set 

at + h 



[0,1] 



fi(x) = 2x — a,; on 



ai, 



, and 



-,bi 



fi(x) = bi on 



As neither nor 1 belongs to K, there is a rational number e > such that 
KC[e,l-e}. Let us define / : [0, 1] -> R by 



We then define t{x) : [0, 1] - 



x on [e, 1 — e], 

2x - (1 - e) on [1 - e, 1] 

e on [0, e] 



by 



t(:c) 



E 2 ^ 



By construction, the function t(x) is computable and non-decreasing, and t(x) > x 
if and only if a; e [0, 1] \ K. As 



i(l)=/(l) = l + e=l + t(0), 



we can take the quotient 



T\ {x) = t(x) mod Z. 



It is easy to see that T\ moves all points towards the set K. More precisely, 
every point iGXis fixed under T\ , and the orbit of every point x <£ K converges 
to inf{y € K n[x, 1]}. Further, for any interval J <g U, all but finitely many T\- 
translates of J are disjoint from J. Hence, no finite invariant measure of T± can 
be supported on J. Thus the support of every Ti-invariant measure is contained in 
K. By Proposition 4.12, no such measure can be computable. 

□ 



We are now equipped to present the counter-example T. We define T 2 : S 1 — > S 1 
by T2(x) — 2a; mod Z, and set 

T = T 1 x T 2 . 

Firstly, by the same reasoning as above, every invariant measure \i of T is supported 
on K x S 1 and hence is not computable by Proposition 4.12. On the other hand, 
T possesses invariant measures of maximal entropy. Indeed, let v be any invariant 
measure of T\ and A the Lebegue measure on S 1 . Setting pi = v x A, we have 



^(T) = /ito P (T 1 ) = log2. 
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5. Harmonic Measure 

5.1. Proof of Theorem C. Let us start with several definitions. 

Definition 5.1. We recall that a compact set K C C which contains at least 
two points is uniformly perfect if the moduli of the ring domains separating K are 
bounded from above. Equivalently, there exists some C > such that for any 
x G K and r > 0, we have 

(B(x,Cr)\B{x,r))r\K = $ =*> K c B(x, r). 

In particular, every connected set is uniformly perfect. 

It is known that: 

Theorem 5.1 (see [17]). The Julia set of a polynomial P of degree d > 2 is a 
uniformly perfect compact set. 

Recall that the logarithmic capacity Cap(-) has been defined in Section 2.3. We 
next define: 

Definition 5.2. Let C C be an open and connected domain and set J = dCl. We 
say that H satisfies the capacity density condition if there exists a constant C > 
such that 

C&p(B(x, r) n <9Q) > Cr for all x G dfl and r < r . (5.1) 

We note: 

Theorem 5.2 (see Theorem 1 in [22]). Condition (5.1) is equivalent to uniform 
perfectness ofdtt. 

The celebrated result of Kakutani [14], gives a connection between Brownian motion 
and the Harmonic measure. 

Theorem 5.3. [14] Let K C C be a compact set with a connected complement f2. 
Fix a point x G Cl and let B t denote a Brownian path started at x. Let the random 
variable T be the first moment when Bt hits d£l, and let uj x = loq,,x denote the 
harmonic measure corresponding to x. Then for any measurable function f on dfl, 

J fduj x =E[f(B T )}. 

In [2] the following computable version of the Dirichlet problem has been proved: 

Theorem 5.4. Let K C C be a compact computable set with a connected comple- 
ment fi. Let x G fl be any point in f2 and let B t denote a Brownian path started at 
x. There is an algorithm A that, given access to K , x, and a precision parameter 
e, outputs an e / '4- approximate sample (3 E from a random variable Bt c where T £ is 
a stopping rule on B t that always satisfies 

e/2< dist(B Tc ,K) < e. 

In other words, we are able to stop the Brownian motion at distance ~ e from the 
boundary. We now formulate the following proposition that is a reformulation of 
Theorem C: 

Proposition 5.5. Let CI C C be the complement of a computable compact set K 
and xq be a point in Q. Suppose CI is connected and satisfies the capacity density 
condition. Then the harmonic measure ujQ tXo is computable with an oracle for xq. 
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Proof. Fix xq 6 O. As before, we denote by (Bt) the Brownian Motion started at 
xq and set 

T = min{t : B t € <9D}. 

We will use Theorem 5.4 together with the capacity density condition to prove 
Proposition 5.5. 

The capacity density condition implies the following (see [12], page 343): 

Proposition 5.6. There exists a constant v = v(C) (with C as in the capacity 
density condition) such that for any r\ > the following holds. Let y £ £1 be a point 
such that dist(y,cj) < T), and let B u be a Brownian Motion started at y. Let 

T y := mm{t : B\ e dfl} 

be the first time B v hits the boundary of CI. Then 

P[\B y Ty -y\> 2n] < v. (5.2) 

In other words, there is at least a constant probability that the first point where 
B y hits the boundary is close to the starting point y. 

Now let / be any function on K satisfying the 1-Lip condition. Our goal is to 
compute 

J fdw=E[f(B T )}. 
within any prescribed precision parameter 5. Note that 

E[f(B T )]=E BT jE[f(B T )\B Te \}]. 
We first claim that we can compute an e such that 

\f 2e (B Te )-E[f(B T )\B Te \}\ < 5/2. (5.3) 

Here T e is given by the any stopping rule as in Proposition 5.4, and /2e(-BT e ) is any 
evaluation of / in a 2e-ncighborhood of Bt c (note that / itself is not defined on 
B Tc i K). 

Let M be a universal bound on the absolute value of /. By (5.2) we can compute 
an e < 5/20 such that if y is e-close to K, the probability that \Bj, — y\ > 5/10 is 
smaller than 5/10M. We split the probabilities into two cases: one where Bt stays 
5/10-close to Bt e , and the complementary case. By (5.3) we have 

\f2e(B Te )-nf(B T )\B Te \]\ = 

\he(B Te ) - E[f(B T )\B Te \, \B T -B Te \< 5/10} \ ■ F[\B T - B Te \ < 5/10] + 
\he{B Te ) - E[f{B T )\B Te \, \B T -B Te \> 5/10] | ■ ¥[\B T -B Te \> 6/10] < 

(5/6) ■ 1 + M ■ (5/10M) < 5/2. 

To complete the proof of the proposition it remains to note that given a (3 e that 
e/4- approximates Bt e as in Theorem 5.4, we can evaluate f2s(Br e ) by evaluating 
/ 3e / 2 (A;) (by evaluating / at any point in a 3e/2-ncighborhood of /3 S ). In this way, 
we obtain 

E^|/ 3£/2 (/3 s )-E[/(B T )]|< 5/2. (5.4) 
Thus, being able to evaluate f^ e /2{Pe) with precision 5/2 suffices. □ 
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5.2. A counter-example. As demonstrated by the following example, even for a 
computable regular domain, the harmonic measure is not necessarily computable. 
Thus the capacity density condition in Theorem C is cruicial. 

For a, b G K, we denote by [a, b] the shortest arc of the unit circle between e 2ma 
and e 2mb . As before, let Ai be the Godcl ordering of algorithms. Define a collection 
of subsets of the circle as follows. If A n halts in I steps, we set j = max(i, 8n) and 
denote by 

L„ := LP n := U^" 2 ^- 1 .^ [2- n + k ■ 2 - 2n -\2~ n + k ■ 2- 2n ^ + exp (-2 2n+2 % 
Otherwise, if A n does not halt, we denote 

L n := L™ := [2~ n - 2" 2 " + 2" 10 ", 2"" + 2" 2 " - 2" 10 "] 
(see Figure 4). 




Figure 4. Left: a set corresponding to an algorithm A n which 
does not halt; right: a set L J n corresponding to A n which halts in 
j steps. 

Let D n denote the disk whose diameter is given by the points cxp(2-7ri(2~" — 
2" 2 ™)) and exp(27h(2-" + 2- 2 ' 1 )). Let 

OO 

A := DU |J D n . 

n=10 

The domain O is obtained by removing the arcs from L n from A. To be precise, set 

OO 

Q:=A\ |J L n . 

n=10 

and 

K :=C\Q. 

We note: 

Proposition 5.7. The compact set K is computable. 
Proof. Note that 

dist ff (^,L^)<12-2- 2 "^. 
We can thus compute L n with an arbitrary precision by emulating A n for sufficiently 
many steps. 

To compute the set K with precision 2~ m , it suffices to compute the first 2 m+A 
sets L n with precision 2~( m+4 \ □ 
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Now let us show that: 
Proposition 5.8. The harmonic measure 

u ■= wn, 

is not computable. 

For a set K C {z : < S < \z\ < r < 1} set 

7(Jf ) :=-logCap(A' ). 
We need an auxilliary lemma: 

Lemma 5.9 (Theorem 5.1.4 in [23]). If K±, . . . , K n are compact subsets of the unit 
disk, then 

1 < 1 1 

7(Jfi U • • • U A„) ~ l{Ki) j(Knj' 

Let S n be the part of the boundary of the disk D n lying outside of D, S n := 
dD n \ D. Harmonic measure is always non-atomic, so if uj is computable, then 
uj(S n ) is also computable. We show: 

Proposition 5.10. If A n does not halt, then Uj(S n ) < 2 _9n+2 . If A n halts, 

Proof. As before, let B t be the Brownian motion started at and let T denote the 
hitting time of dfl, 

T := inf{i : B t G <9ft}. 
Let us recall that for E C 5il we have 

uj{E) = P[B T G A]. 

Assume now that A n docs not halt. Let us introduce a new domain tt' := 
C \ [2"" - 2- 2 ™,2~™ + 2~ 2 ™] and T" be the corresponding hitting time. Observe 
that if By G S« then 

£?jv G A„ := [2"" - 2~ 2? \ 2- n - 2~ 2n + 2~ 10n ] U [2"" + 2~ 2rl - 2" 10 ", 2"" + 2~ 2 "]. 
Thus 

un,o(Sn) = P[St G £„] < P[fl T , G A„] = l* () o(An). 
The desired estimate is now obtained by mapping (f2',0) conformally to (D, 0). 

Assume that A n halts in j steps. To bound co(S n ) in this case, we will use the 
following estimate on harmonic measure ([12], Equation (III. 9. 2)): 

Let A C {z : < 5 < \z\ < r < 1}. Then 



log (V) 
>y{K ) + log(l - r 2 ) 



w D \K ,o(A ) < ■ ^ (5.5) 

Let T" denote the hitting time of <9B by B t , and let 

M n := {z G L~ : dist(z,0ft) > 2- 2 "-^ 4 } 

be the part of the arc lying relatively far away from the boundary. 

Conformally mapping D n to the unit disk centered at z € M n and using the 
estimate (5.5) and Lemma 5.9, we obtain that for z G M n we have 

u Dn \ Ln A L n) < 2" J+3 < 1/8 (5.6) 
We will also need T\ > T" - the first hitting time of dD n after hitting dD, and 
T2 < Ti - the first hitting time of dD n U L n after hitting dD. 
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Note now that 

F[B T „ G M„] = lcngth(M„)/27r > (5.7) 
Let us note that by symmetry and estimate (5.6), we have 

¥[B T g 5 n | B T „ g M„] > P[B Tl g s„ | s T » e M n ]-V[B T2 g L„ | 5 t - e M n ] - 

i - ¥[B T , G L n | i?T" G M n ] > - - max Wi, B \£ ni ,(Z n ) > 1/4. (5.8) 

Z Z z£M n 

The desired lower estimate on cj is now obtained by combining (5.7) and (5.8). 

□ 

We now conclude the proof of Proposition 5.8: 

Proof of Proposition 5.8. Assume the contrary, that is, suppose that ui is com- 
putable. For every n G N let {Uj l (z)}°° =1 be a sequence of functions given by: 

C 1 if dist(z,S„) < 2~ j ; 
Uf{z) = I if dist(z, S„) > 2 • ; 

[ 1 - 2J {d - 2~ j ) if d = dist(z, S n ) G [2~ j , 2 ■ 2~ j ] 

We have: 

(a) the functions U™(z) are computable uniformly in n and j. 
Since u> is non-atomic, 

(b) for a fixed n we have 

m (z)du > w(5«) and / U" {z)duj — w(5«). 

Similarly, we can costruct a sequence of functions L™(z) such that 

(c) the functions £™(z) are computable uniformly in n and j; 

(d) for a fixed n we have 

L r Hz)du < u){S n ) and / L n Az)du — -> w(5„). 

We leave the details of the second construction to the reader. 

By part (3) of Proposition 3.10, properties (a) and (c) imply that the integrals 

UJ(z)du and J L](z)duj 

are uniformly computable. Consider and algorithm .Ahait which upon inputting a 
natural number n does the following: 

(1) 3 ■= i; 

(2) evaluate Uj,lj such that 
U?(z)du\ < 2- 20n and \lj - J L]{z)dcj\ < 2" 20 "; 

(3) if Uj < 2" 9n + 2 + 2- 19n then output and halt; 

(4) if l 3 > 2- 2n - 3 - 2- 19n then output 1 and halt; 

(5) j :=j + l and go to (2). 
By Proposition 5.10 and properties (b) and (d), we have the following: 

• if An halts then .Ahait outputs 1 and halts, and 

• if An does not halt then .Ahait outputs and halts. 
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Thus v4hait is an algorithm solving the Halting Problem, which contradicts the 
algorithmic unsolvability of the Halting Problem. □ 
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